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We investigate the influence of disorders in terminal ethylene groups of BEDT-TTF 
molecules (ethylene-disorders) on the vortex pinning of the organic superconductor k-(BEDT- 
TTF)2Cu[N(CN)2]Br. Magnetization measurements are performed under different cooling- 
processes. The second peak in the magnetization hysteresis curve is observed for all samples studied, 
and the hysteresis width of the magnetization becomes narrower by cooling faster. In contradiction 
to the simple pinning effect of disorder, this result shows the suppression of the vortex pinning 
force by introducing more ethylene-disorders. The ethylene-disorder domain model is proposed for 
explaining the observed result. In the case of the system containing a moderate number of the 
ethylene-disorders, the disordered molecules form a domain structure and it works as an effective 
pinning site. On the contrary, an excess number of the ethylene-disorders may weaken the effect of 
the domain structure, which results in the less effective pinning force on the vortices. 



I. INTRODUCTION 

BEDT-TTF based organic superconductors, where 
BEDT-TTF is bis(ethylenedithio)tetrathiafulvalene, 
have the layered structure and the two-dimensional 
electronic properties^ As a feature of molecular-based 
compounds, the internal degree of freedom in the 
BEDT-TTF molecule can give a positional disorder. 

The terminal ethylene groups of the molecule can be 
thermally excited between two conformations, "eclipsed" 
and "staggered" (Fig. [Q. In the present compound, k- 
(BEDT-TTF) 2 Cu[N(CN)] 2 Br, one of the ethylene group 
in a molecule is known to be disordered and the other 
ordered at room temperature^ where the eclipsed con- 
formation is stable. Such a thermal excitation is ob- 
served in the narrowing effect of the 13 C-NMR line width 
above ~150 K£ The conformational disorder (ethylene- 
disorder) is frozen by quenching, and the introduction of 
the ethylene-disorders slightly affect the electronic state. 
The resistivity below about 80 K depends on cooling- 
rate^ and in addition the superconducting transition 
temperature T c is also reduced by cooling faster^ The 
heat capacity measurement suggested that a glass tran- 
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FIG. 1: Schematic view of conformational disorder in BEDT- 
TTF. 



sition is caused by freezing of the motion in the ethylene 
group around 100 K& Tanatar et al. proposed the T-P 
phase diagram on ordering of the ethylene groups^ Stal- 
cup et al. reported that a time-dependent resistivity is 
correlated with T c and the amplitude of the Shubnikov- 
de Haas oscillation in terms of the ethylene-disorders >2> 
However, none of the reports have clarified the spatial 
distribution of the ethylene-disorders at all. 

On the superconducting mixed states of the K-phase 
BEDT-TTF organic superconductors, following features 
have been observed; the wide reversible magnetization 
regime^ the second magnetization peaksfi2«i*i2iia4ii and 
the first-order phase transition between the vortex solid 
and liquid states>i^*iSiiLi^ Taniguchi and Kanoda re- 
ported the cooling-rate dependence of the magnetization 
curves^ where the hysteretic magnetization width be- 
comes narrow by cooling fast. One may consider that 
cooling fast introduces a much larger number of the 
ethylene-disorders, and leads to much stronger vortex 
pinning on the basis of a concept of the simple pinning 
effect by disorder. Nevertheless, this expectation does 
not agree with their experiments. The nature of this 
controversial result is unclear. 

In this paper, we report detailed measurements of the 
hysteretic magnetizations under the control of cooling- 
process. We clarify the effect of the ethylene-disorders on 
the vortex pinning, and discuss the spatial distribution of 
the ethylene-disorders for explaining the observed results. 
We propose that domain structures originated from the 
ethylene-disorders work as a pinning site. 

II. EXPERIMENTALS 

Single crystals were synthesized by means of electro- 
crystallization. All the measurements of the magneti- 
zation were performed using a SQUID magnetometer 
(Quantum Design, MPMS-5) for a single crystal with 
the dimension of 1.8 x 1.6 x 0.3mm 3 . No grease was 
used to avoid strain to the sample. The magnetic-field 
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FIG. 2: Temperature dependence of the static magnetic sus- 
ceptibility at 5 Oe. The closed and open symbols indicate 
the shielding curve measured under a zero-field-cooling con- 
dition (ZFC) and the Meissner curve under a field-cooling 
(FC), respectively, with overlaid for the slow-cooled (circles), 
rapid-cooled (triangles), and quenched processes (squares). 
The broken line shows the full-Meissner volume. 

was applied perpendicular to the conduction plane. The 
cooling-process was controlled from room temperature to 
15 K as follows; (1) 0.17 K/min with 30 h anneal at 70 
K (slow-cooled), (2) 15 K/min (rapid-cooled), and (3) 
100 K/min (quenched). After all the sequential mea- 
surements of the quenched sample, the data taken in the 
slow-cooled process again were completely reproduced, 
suggesting no deterioration of the sample by quenching, 
such as micro cracks. We carried out the similar mea- 
surements for other two samples, which indicate the same 
results with the present report. 

III. RESULTS 

Figure [21 shows the temperature dependence of the 
static magnetic susceptibility at 5 Oe. After subtracting 
the contribution of the core diamagnetization (—4.8 x 
10~ 4 emu/mol), the demagnetization factor was cor- 
rected using an ellipsoidal approximation. T c 's defined 
as an intercept of the extrapolated lines of the normal 
and reversible susceptibilities are 11.9, 11.5, and 11.0 K 
in the slow-cooled, rapid-cooled, and quenched process, 
respectively (short arrows in Fig. The suppression 
of T c by cooling faster is quantitatively consistent with 
the previous report^iS The magnetizations under a zero- 
field-cooled (ZFC) condition (shielding curves) show al- 
most the full Meissner volume and the negligible cooling- 
rate dependence within the accuracy of our measure- 
ment. This result demonstrates that the superconducting 
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FIG. 3: Magnetic-field dependence of (a) the magnetization 
at 5.5 K, and the magnetization hysteresis width for cooling- 
rates of (b) slow-cooled, (c) rapid-cooled, and (d) quenched 
condition, respectively. 



volume fraction is almost unchanged by cooling-process, 
whereas Pinteric et al. reported a small suppression of 
the superconducting volume fraction by faster-cooling. 19 
On the other hand, the magnetic susceptibilities under a 
field-cooled (FC) condition (Meissner curves) are reduced 
in comparison with those under the ZFC condition (open 
symbols in Fig.EJl. The volume fractions are slightly sup- 
pressed with increasing cooling-rate. At first sight, the 
suppression might be attributed as reinforcement of vor- 
tex pinning by cooling fast, because smaller exclusion of 
magnetic-field generally reflects stronger pinning. How- 
ever, as shown later, the experimental results in higher 
magnetic-fields are inconsistent with this explanation. 
We so suggest the other scenario as follows. It should be 
noted that the magnetic-field of 5 Oe for the susceptibil- 
ity measurements is not sufficiently smaller than H c i(0) 
= 19 G. 20 -ff c i(0) seems to decrease with T c by cooling 
faster, and therefore the effective field H/H c i(Q) is in- 
creased. This makes the flux exclusion more incomplete, 
leading to the further suppression of the FC supercon- 
ducting volume by cooling faster. 

Figure^a) shows the magnetic-field dependence of the 
magnetization at 5.5 K. A typical second peak is found 
around 500 Oe. As displayed in the magnetization hys- 
teresis width AM (Fig.|3fb)-(d)), the second peak is ob- 
served below about 8 K for all the cooling- process. The 
peak height seems to decrease with increasing cooling- 
rate as reported previously. 12 Strictly speaking, to depict 
the cooling-rate dependence of AM, a comparison of the 
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measured AM's should be done at a fixed temperature t 
(= T/T c ), since cooling-rate changes T c , 

Figure 0] shows the magnetic-field dependence of the 
critical current density J c at t — 0.4f . The estimation of 
J c from AM was carried out through a standard Bean's 
modelSi for a square-shape sample 22 . The data at 4.5 K 
are used for the quenched sample, and the correspond- 
ing data at 4.7 K (rapid-cooled) and 4.9 K (slow-cooled) 
are calculated from a linear interpolation with the data 
at both 4.5 and 5.0 K, respectively. Even after the sub- 
traction of the effect on the reduced T c , the critical cur- 
rent density clearly decreases with increasing cooling- 
rate. This indicates the suppression of the effective vor- 
tex pinning by cooling fast. Because rapid cooling of 
samples is considered to induce frozen ethylene-disorders, 
one may expect the increase of vortex pinning. However, 
in the present vortex system, the effective pinning force 
becomes weak with increasing ethylene-disorder. This 
result is apparently inconsistent with a concept of the 
simple pinning driven by disorder. 

Here we characterize the magnetic-field of the second 
peak. Contrary to the peak effect of Bi2Sr2CaCu20s 
(Bi2212)f2£ the peak structure in the present supercon- 
ductor is broad. This makes it difficult to determine a 
characteristic magnetic-field of the second peak without 
any ambiguity. We here define the peak position using 
two magnetic-fields: one is the magnetic-field of the mini- 
mum point between the central and second peaks (H m - m ), 
and the other is the maximum point for the peak (H p ). 
Each magnetic-field is plotted in the T-H phase diagrams 
with overlaid for each cooling-process, as shown in Fig. 
|SJ Additionally plotted irreversibility fields (-ffi rr , open 



FIG. 5: Characteristic magnetic- fields (closed symbols) as a 
function of temperature with overlaid for the slow-cooled (cir- 
cles), rapid-cooled (triangles), and quenched (squares): (a) 
H p and (b) Hmin- H[ TI (open symbols) is also plotted. 



symbols) are defined using a criterion of J c = 2.0 A/cm 2 
(corresponding to AM = 1.0 X 10~ 5 emu). As a function 
of temperature, H p increases with reducing temperature 
for all the cooling-process (Fig. 0a)). On the contrary, 
-ffmin shows the different behavior as displayed in Fig. 
EJb), where the H mm lines are piled around 200 Oe. 



IV. DISCUSSION 

We first discuss the origin of the second peak ef- 
fect. The peak effect of the present compound has been 
attributed to the dimensional crossover of the vortex 
systemiiSiiiiiS A vortex line changes to two-dimensional 
(2D) pancake vortices with the decrease of the inter- 
layer Josephson coupling above the crossover field B cx sa 
4>o/{jd) 2 , where <j>o is the magnetic flux quantum, 7 the 
anisotropy factor, and d the interlayer distance. With 
the values of 7 = 130 - 18Q2i2£ and d = 15 A, one 
obtains B cr k 280 - 540 G. The estimated B cr roughly 
corresponds to the second peak. It is almost the same 
as -ffmin and H p , and thus the dimensional crossover sce- 
nario as the second peak origin is applicable to this case. 
Although the broad shape of the second peak seems to 
be inconsistent with this interpretation, the broadness 
may come from the averaging of the peak structure over 
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the bulk sample in which the local magnetic-field has 
a distribution. A more detailed experimental investiga- 
tion would be needed to check this point. B cr is gen- 
erally independent of temperature, and moreover in the 
present system, cooling-rate independence is expected as 
well, because the anisotropy factor 7 is likely to be un- 
changed by cooling-rate. Actually, even in the deuterated 
BEDT-TTF based-salt, which has dramatically different 
electronic properties from the present (non-deuterated) 
saltr& the same value of 7 has been reported^ While H p 
decreases with increasing both temperature and cooling- 
rate, -ffmin is almost independent of both of them. The 
behavior of H p cannot be explained by a simple dimen- 
sional crossover model. On the contrary, i? m ; n seems to 
be appropriate as B cr . If B CI ~ i? m in, the second peak is 
qualitatively explained as follows: In H >^ H m i n , two- 
dimensional pancakes appear, and the pinning force for 
the pancake vortices is effectively optimized around H p . 

Besides dimensional crossover, several models have 
been proposed as the origin of the second peak effect: 
(1) matching effects^ (2) the magnetic-field-dependent 
relaxation-rate (collective pinning model) (3) the field- 
induced order-disorder transition (Bragg-glass to vortex- 
glass transition), etc. First, in (1) matching effects, the 
peak structure appears around a magnetic-field in which 
the vortex lattice structure is commensurate with the 
defect structure. In this model, when the number of 
defects increases, the magnetic-field of the peak maxi- 
mum increases. However in the present case, by cool- 
ing faster which may be giving much larger number of 
defects, H p shifts toward lower magnetic-fields. Thus 
the experimental results are inconsistent with this model. 
Next, the second model (2) is also excluded by the fol- 
lowing reason; in this model, a collective pinning of vor- 
tices appears above the peak maximum field, resulting in 
a slower relaxation-rate than that below the magnetic- 
field where a single-vortex pinning mechanism governs 
the system. According to a relaxation measurement in 
the present compound^, the relaxation-rate at the peak 
maximum field is faster than that at the peak minimum 
one, which contradicts with this model. The third model 
(3) is relatively controversial. In this model, with increas- 
ing magnetic-field, a first-order phase transition happens 
from a vortex lattice or Bragg glass state (ordered-solid 
phase) to a vortex-glass state (disordered-solid phase), 
lacking translational order. This transition is originated 
from the larger pinning energy gain of the glass phase 
than that of the lattice one, which can reasonably ex- 
plain the peak effect. In high-T c superconductors^^!, 
the order-disorder transition has been observed as a steep 
increase or a step of magnetization, but in the present 
organic superconductor, there has been no decisive evi- 
dence in favor of this model. To confirm this model, it 
is needed to investigate, for example, the magnetization 
under a rapidly fluctuating magnetic-field, as has been 
measured in Bi2212i2i 

While the dimensional crossover explains the second 
peak effect, the problem of the suppressed J c by cooling 



fast is left unsettled. Note that the suppression is remark- 
ably shown in H > H m i n (Fig. 0J , in which the vortex 
system consists of pancakes. This implies the weakness 
of pinning, which will support the discussion on domain 
structures of the ethylene-disorders as a weak pinning 
origin, presented later. In the following, to understand 
the observed result we solve the puzzling origin of the 
vortex pinning. 

As a pinning site, intrinsic defects such as molecular 
defects, dislocations, or coexist of another crystal phases 
are considerable. Actually a coexisting a-phase in the re- 
type crystal has been observed using scanning tunneling 
microscopy (STM)i22i However, these intrinsic defects, if 
any, are independent of cooling-rate. As already men- 
tioned, our experimental results indicate that the effec- 
tive pinning force is suppressed with cooling faster, that 
is, with increasing ethylene-disorder. Here it should be 
noted that the size of a BEDT-TTF molecule (~ 3 A 
projected on the ac-plane) is much smaller than that of 
a vortex core (2£|| (0) ~ 50 A, where £|| (0) is the in-plane 
coherence Iength2i21). Thus it is difficult for an individ- 
ual ethylene-disorder to behave itself as a pinning site 
because of smallness of the molecular size. 

We next describe the numbers of the ethylene- 
disorders. In the present case, the quenched condition 
will give the largest number of the ethylene-disorders, be- 
cause the conformational disorders at room temperature 
are frozen by quenching. The fraction of the ethylene- 
disorders in the quenched condition is suggested to be 
about 5-20%&^ Although there are no corresponding 
reports after slow-cooling, we suppose that a slight num- 
ber of the ethylene-disorders exist even in our slow-cooled 
condition. Tokumoto et al£ reported that T c monoton- 
ically decreases from 12.2 K to 11.7 K as the cooling- 
rate increases from 0.02 K/min up to 10 K/min. T c of 
11.9 K at the rate ~ 0.2 K/min in their report shows a 
good agreement with our slow-cooled case. Thus T c of 
the present sample can furthermore increase by cooling 
more slowly. In other words, this result suggests that 
there exist some ethylene-disorders even in the present 
slow-cooled condition. 

Next, we discuss the spatial distribution of the 
ethylene-disorders. At high temperatures, thermal ex- 
citation will cause a random mixing of the staggered and 
eclipsed molecules and/or these micro domains. The 
quenched process freezes a higher temperature glassy 
state with a random distribution of the disordered 
molecules. As the temperature is slowly lowered, the 
stable conformation (eclipsed) gradually becomes domi- 
nant. The remnant ethylene-disorders weaken the super- 
conducting state with reducing T c . Although the mecha- 
nism in the relaxation process of the ethylene-disorders is 
unclear, several relaxed states after slow-cooling are ex- 
pected. The possible relaxed states are classified into fol- 
lowing three states as shown in Fig.[S] As an example, we 
here assume 5% concentration of the ethylene-disorders, 
which will be the maximum in the slow-cooled condition. 6 
Each panel shows the molecular arrangement projected 
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FIG. 6: Schematic views of relaxed states in the slow-cooled 
condition with 5% ethylene-disorders. Each panel shows the 
molecular arrangement projected on the ac-plane. Three pos- 
sible distributions of the ethylene-disorders are described: (a) 
a random distributed state, (b) a state with ordered domains 
originated from vortex core-scale inhomogeneity, and (c) a 
state consisting of ordered and dense-disordered domains. 



on the ac-plane over ~100 A, containing about 640 
BEDT-TTF and 32 disordered molecules. The BEDT- 
TTF molecules are described as small circles with diam- 
eter of about 3 A, consisting of ordered (open symbols) 
and disordered (closed symbols) molecules. Figure EJa) 
shows a state with ethylene-disorders distributed ran- 
domly, which is similar with the quenched one except for 
the small number of the ethylene-disorders. The second 
possible state shown in Fig. Efb) features the inhomo- 
geneous distribution of ethylene-disorders. In this case, 
ordered domains are formed by the increase of the stable 
eclipsed molecules, while remnant disorder-including ar- 
eas are coexisting. A boundary structure can be formed 
between the ordered domain and the disorder-including 
area, whereas the spatial structure of the domain bound- 
aries is unidentified. The third state is containing do- 
mains of dense-disordered molecules (Fig.Efc)). The ex- 
istence of domain structures consisting of only staggered 
or eclipsed molecules were suggested by the study on the 
relaxation of the resistivityiM Here it should be noted 
that the interlayer correlation of the disorder's distribu- 
tion is now out of question, since we now treat the pinning 
of the two-dimensional pancake vortices. Therefore our 
description of the ethylene-disorder's distribution is also 
two-dimensional one, but this does not rule out the ex- 
istence of the interlayer correlation of the domain struc- 
ture. 

Taking into account the above information on the 



ethylene-disorders, we propose the domain structures of 
the relaxed state shown in Fig. EJb) as a vortex pin- 
ning origin. This relaxed state contains the characteris- 
tic boundary structure between the ordered domain and 
the disorder-including area. The size of such a domain 
structure is much larger than at least the in-plane coher- 
ence length ~ £|| (the vortex-core-scale). In the slow- 
cooled condition, the widely-ordered domains dividing 
the disorder-including areas partly are formed, result- 
ing in the vortex-core-scale inhomogeneity. Therefore, if 
these domain structures work as vortex pinning sites, the 
slow-cooled process gives much stronger effective pinning 
than the quenched one. 

On the other hand, in the relaxed state shown in 
Fig. EJa), many individual ethylene-disorders are dis- 
tributed randomly, giving a quasi-homogeneous state on 
the vortex-core-scale. Therefore, the ethylene-disorders 
do not work as vortex pinning sites, whereas they give rise 
to many point defects. The quenched condition is similar 
with this state, while the concentration of the ethylene- 
disorders is much larger. Such an excess number of the 
ethylene-disorders may weaken the effect of the domain 
structure, resulting in the less effective pinning force on 
the vortices. This is consistent with our experimental 
results. 

One may expect that the dense-disordered domain 
shown in Fig.[SJc) gives a strong pinning origin. However, 
the size of the domain, in spite of containing 5% ethylene- 
disorders, is still smaller than the vortex core. Of course 
one can describe a much larger dense-disordered domain 
structure in 5% ethylene-disorders than that of Fig.|Bfc), 
but then the total number of the domains over the sam- 
ple size becomes quite a few. So the contribution to the 
vortex pinning by such a dense-disordered domain cannot 
be so strong, if any. 

An experimental support for the existence of a spa- 
tial inhomogeneity was reported by Tanatar et al2^, who 
observed a magnetic viscosity phenomena and proposed 
a magnetic domain structure formation in the antiferro- 
magnetically ordered state in the normal state. Unfortu- 
nately, it is unclear whether their magnetic domain struc- 
ture is directly related to the ethylene-disorders. While 
this may be a possible solution to the present problem, 
further experimental studies will be needed to check this 
point. 

A clue might lie in a direct observation of the domain 
structures with local probe, such as STM, which is plan- 
ning for the near future. 



V. CONCLUSIONS 

In conclusion, we reported the hysteretic magnetiza- 
tion measurements to investigate the disorder effect on 
the vortex pinning in the cooling-rate of 0.17, 10, and 
~ 100 K/min. The negligible cooling-rate dependence 
of the almost full volume fraction in the superconduct- 
ing state demonstrates that the superconducting volume 
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is unchanged by cooling-process. The second peak is 
observed below about 8 K for all the cooling-process. 
We adopt the dimensional crossover scenario as the ori- 
gin of the second peak. The vortex pinning in the 
present system is suggested to be closely related to the 
ethylene-disorders. However, it is difficult for an individ- 
ual ethylene-disorder to behave itself as a pinning site 
because of smallness of the size. Actually by cooling 
faster the critical current density decreases, and so the 
effective pinning force is apparently suppressed. Thus a 
simple pinning by the introduced ethylene-disorders does 
not occur. 

To explain the puzzling pinning origin, we proposed the 



ethylene-disorder domain model, where domain struc- 
tures form a vortex-core-scale inhomogeneous distribu- 
tion of the ethylene-disorders working as vortex pinning 
sites (shown in Fig. E^b)). The size of the domain struc- 
ture will be much larger than at least the in-plain co- 
herence length. On the contrary, an excess number of 
disorders introduced by quenching may weaken the ef- 
fect of the domain structure, which results in the less 
effective pinning force. 

This research was partially supported by the Ministry 
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for Encouragement of Young Scientists. 
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